Prom Bubble to Skyrmion: Dynamic Transformation Mediated by a Strong Magnetic 

Tip 



m 
O 
o 

(N 



(N 



A.S. Kirakosyan 1 

Department of Physics, Texas A&M University, College Station, Texas 77843-4242 

V. L .Pokrovsky 1 » 2 

1 Department of Physics, Texas A&M University, College Station, Texas 77843-4242 
2 Landau Institute of Theoretical Physics, Moscow, Russia 
(Dated: February 2, 2008) 

Skyrmions in thin metallic ferromagnetic films are stable due to competition between the RKKY 
interaction and uniaxial magnetic anisotropy. We mimic the RKKY interaction by the next-nearest- 
neighbors ferromagnetic and antiferromagnetic exchange interactions. We demonstrate analytically 
and numerically dissipative transformation of a bubble created by a strong cylindrical magnetic tip 
into a stable Skyrmion. 
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1. INTRODUCTION 

Recently a significant progress has been achieved in the 
nanofabrication and design of new technologies for very 
high density storage on hard disks (densities above 500 
Gbit/in 2 )i and magnetoresistive random access mem- 
ory (MRAM) . For the further miniaturization of the data 
storage and the magnetic memory it is important to de- 
crease the size of nanoelements. Magnetic defects in such 
nanoelements storing bits of information can have very 
small size. The properties of Skyrmion in thin metallic 
magnetic films (the small size about 1 nm and stability) 
may allow to use their two-dimensional (2D) arrays in 
magnetic storage devices. To the best of our knowledge 
a magnetic memory based on Skyrmions was not devel- 
oped yet. There are no analogs of storage of so extremely 
high density. 

Skyrmions are strongly localized scale invariant 
analogs of magnetic bubble domains (MBD) 2 . The axi- 
ally symmetric Skyrmion solution for boson fields mini- 
mizing the Hamiltonian and describing the structure of 
nuclear matter was first found by Skyrme 3 .. Belavin and 
Polyakov based on nonlinear cr-model have established 
that Skyrmion is a topologically nontrivial metastable 
minimum of energy - *. Since then the properties of 
Skyrmions in two-dimensional magnets and thin mag- 
netic films have been widely investigated in last few 
decadea 5 . 6 . 7 . 8 . 9 . 

Slow relaxation spin dynamics of Skyrmion was stud- 
ied by Abanov and Pokrovskji They calculated the 
radius of the Skyrmion as a function of time. They 
have found that the interaction of the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) type together with the uniaxial 
magnetic anisotropy can stabilize the Skyrmion. 

In quasi-2D antiferromagnets the Skyrmions were in- 
directly observed by F. Waldner by the ESR 1 — and the 
heat capacity measurements^ in a good agreement with 



the Belavin and Polyakov (BP) theory. The appearance 
of Skyrmions in the QHE was predicted by Sondhi et 
al4£. In the regime of fractional QHE Skyrmions were 
discovered experimentally in GaAs quantum wells by 
means of optically pumped NMR and by measurements 
of the Knight shift, which is proportional to the elec- 
tron magnetization^. For integer QHE Skyrmions were 
detected in thermal magnetotransport measurements 14 
and in optical absorption 1 ^. 

Recently another topological excitations, magnetic 
vortices were observed and studied experimentally in 
permalloy and Co nanodisks by means of the Lorentz 
microscopy and Magnetic Force Microscopy (MFM) 16 . 
The evolution and stability of such magnetic vortices in 
a small isotropic cylindrical magnetic dot was analyzed 
theoretically ini£. 

The equilibrium size of Skyrmion in 2D thin magnetic 
films is extremely small. This is one of serious obsta- 
cles for its experimental observation. For soft magnets 
the size of vortices is defined by competition between 
the exchange and magnetic dipole interaction. In uni- 
axial magnets, strong uniaxial anisotropy together with 
the fourth order exchange interaction leads to small size 
of Skyrmion. Recent great progress in nanotechnology 
gives a chance to detect experimentally Skyrmions in thin 
magnetic films. The experimental methods employed for 
studying Skyrmions in QHE systems and vortices in soft 
magnetic dots can be also applied to 2D uniaxial mag- 
nets. So far the only technique potentially capable of 
atomic resolution is the spin-polarized scanning tunnel- 
ing microscopy (SPSTM). 

Due to stability of the Skyrmion in an uniaxial thin 
metallic magnetic film it is possible to create it by a 
strong magnetic tip. In this paper we study the prop- 
erties of the MBD created by the tip as well as its re- 
laxational dynamics when tip moves away uniformly and 
vertically from the film. When the distance between the 
film and the tip exceeds a critical value, it is not able to 
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confine the MBD anymore and the latter collapses and 
gradually transforms into a stable Skyrmion. For numer- 
ical simulations we use a stabilization mechanism pro- 
vided by the interplay between next nearest neighbor an- 
tiferromagnetic interaction and the magnetic anisotropy. 
We also show that the demagnetization energy of the 
Skyrmion is equivalent to the renormalization of the sin- 
gle ion anisotropy constant. We find that demagnetizing 
field produced by a single Skyrmion is strongly localized 
around Skyrmion center. 

The plan of this paper is as follows. In Section [21 we 
introduce a theoretical model and equations of motion. 
We calculate the demagnetization energy of the Skyrmion 
in Section H Section 01 is devoted to the magnetic field 
generated by a cylindrical tip. In Section we study 
static properties of the MBD and the Skyrmion solv- 
ing numerically Landau-Lifshitz equation with the fitting 
method. In Section [f)| we study analytically and numer- 
ically a discrete spin model with the dipolar interaction 
and find stability condition for the static Skyrmion. Re- 
laxational dynamics of the transformation of a MBD into 
a Skyrmion in the presence of an uniformly moving tip is 
considered in Sectional Our conclusions and discussions 
are presented in the last section [SJ 



2. MODEL 

We consider the two-dimensional Heisenberg model 
with the exchange interaction up to the fourth order 
derivatives and uniaxial anisotropy in continuous approx- 
imation. Its Hamiltonian has a following form 

H = h J W(M)d 2 x (1) 
where W(M) = W e (M) + W a (M) + W Z (M) and 

r!VM) 4 ) (2) 
(3) 



w e (M) = - ^(vm) 2 , 

W a (M) = (W-M 2 ), 



Wz{M) = H Z M - (^H(r,d(t)) + iH D ^ M, (4) 

Here the energy density iy(M) consists of the exchange 
energy W e (M.), the uniaxial anisotropy energy W a (M.); 
Wz(M) contains the interaction energy of cylindrical tip 
with the magnetic film and the demagnetization energy. 
In our notations a is the exchange interaction constant, 
[3 is the easy- axis anisotropy constant, n is a positive 
dimensionless constant, characterizing the fourth order 
exchange interaction (we consider just this term because 
of its simplicity), Mo is the spontaneous magnetization 
of the film, H(r, d) is the magnetic field created by the 
magnetic tip (cylinder) at the surface of the film as a 
function of the polar radius r and the distance d between 
the cylinder and the plane of the film; h is the film thick- 
ness. In a thin magnetic film the demagnetization energy 



of a localized magnetic excitation with the radius R e at 
R e <C h simply renormalizes the anisotropy constant (3 
(/3 — >■ ft — 47r) . In the opposite case R e > h the demag- 
netization energy is more complicated. Below we calcu- 
late it for the Skyrmion. For iron the exchange length 
lex = 1.5nm. Therefore, an ultrathin iron film consists 
of few monolayers. 

The dynamics of a classical 2d ferromagnet is governed 
by the Landau-Lifshitz equation (LLE) for the magneti- 
zation vector M 



MM) 

rM(r,() x M(r,f) 
Mo 



7M(r,t) x 



5H[M] 
SM(r,t) 



(5) 



here 7 = g/is/ft; g and [Ib are the gyromagnetic ratio 
and the Bohr magneton, respectively, and T is the di- 
mensionless Gilbert damping constant. The LLE equa- 
tion JSJ in angular variables 9 and <p can be derived from 
the Lagrange-Hamilton variational principle with the fol- 
lowing Lagrangian density C{6, 4>) 



C(9,, 



Mi 



(1- cos0) -Z-W(0,<l>) (6) 
7 at 



In the static case the LLE Eq. (J5J read: 



Mn 2 J^=M(r) (M(r) ™ 



<JM(r) 



<5M(r) 



(7) 



The right-hand side of Eq. J7J originates from a con- 
straint: M 2 = Mq at any point r. 

The classical ground state is double degenerate due to 
uniaxial magnetic symmetry. A topologically nontriv- 
ial minimum of the Hamiltonian Q corresponds to the 
Skyrmion solution 



Ml = M, 



■ 



i? 2 + r 2 
r 2 -Rl 
■ 2 +R 2 ' 



(8) 



The Skyrmion is characterized by its radius R s and an 
angle ifi of the global spin rotation; r and tp are polar co- 
ordinates. For isotropic magnets the energy of Skyrmion 
does not depend neither on R s nor on ip. In real systems 
the degeneracy with respect to R s is lifted by uniaxial 
magnetic anisotropy and by the fourth order exchange 
interaction. 

Another static solution of the LLE is the well known 
MBD, which was the subject of intensive study during 
last few decades. Here we consider an MBD with the 
topological charge one. The structure of such a soliton 
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FIG. 1: Demagnetizing field of the Skyrmion 1131 versus po- 
lar radius normalized to the film thickness r/h. Solid and 
dashed lines correspond to Hj 3 /(4-kMo) and / '(Air Mo), re- 
spectively. The radius of Skyrmion R s = 2h. 



can be approximately described by the "domain- wall" 
ansatz: 



M b x 


= Mo 


sech— 


K 


= M 


sech— 


Ml 


= M 


tanh 



cos(y + tp), 



h 

r-R b 
lo 



sm(tp + ip), 



(9) 



where Rb is the MBD radius. Since we consider an ultra- 
thin magnetic film, the MBD exists only in the presence 
of an external magnetic field and collapses without this 
field. 



3. DEMAGNETIZATION FIELDS 

Let us now return to the demagnetization energy. The 
magnetostatic problem for the non-topological MBD was 
solved first by Thiel o 18 ! 19 . He expressed Hd in terms of 
elliptic integrals. The Thiele's solution for demagnetiz- 
ing field x~ 1 (d£n /dx) can be approximated with surpris- 
ingly good accuracy 



x- 1 (d£ D /dx) « -(1 + &E/2)- 



(10) 



where dimensionless radius x — Rb/h and dimension- 
less energy £d = (16ir 2 M{j h 3 )~ 1 H.D are introduced. For 
magnetostatic interaction the topological charge does not 
play any role. Hence the Thiele's solution persists when 
dipolar forces are taken in account. 

For ultrathin magnetic films consisting of few layers the 
ratio of the exchange length to the thickness is l ex /h > 
4/3. Therefore, the MBD solution is unstable without 
external magnetic field. 



The reading of the information stored by Skyrmions 
can be performed by the measurement of the Skyrmion 
magnetic field. Below we calculate this field and its en- 
ergy. The Fourier transform of the Skyrmion magnetic 
field has the following form 



h D (k) 



-47rk 



(k-m(k)) 

k 1 



(ii) 



The magnetostatic energy in terms of the magnetization 
Fourier-components reads: 



En = 



1 



(27T)* 



k • m(k) 

k 1 



-dk 



(12) 



To calculate the magnetic field we substitute the Fourier 
transform of the Skyrmion magnetization JSJ into {TJ| 
and perform the inverse Fourier transformation. The re- 
sult is as follows: 



iff = 4ttM R 2 s 



*? = 
where 

F z = 

F r = 



4ttM Rt 



F z (q, z, h) K (qR s ) J (qr) qdq 
-4nM &[h/2- \z\] 
F r (q, Z, h) K (qR s ) J x (qr) qdq, (13) 



sign 



sign 



+ z 



,-«*+* 



(14) 



The origin of the frame of reference coincides with the 
Skyrmion center. The Skyrmion demagnetizing field for 
R s /h = 2 is shown in Fig. ^ The demagnetization 
energy of the Skyrmion l|12l) in units Eq has the following 
form: 



Er 



E D -jE a 



(15) 



where the last term is a part of magnetostatic energy 
renormalizing anisotropy constant: j3 — > j3 — 47r 6 and £p 
has the form 



f 1 
t n 



np° 



18 

,-,2 



(Stt-AE (m)-2K (to)) 
p F q ({1,1, 3}, {5/2, 5/2}, to) 



+ "- p F q ({1,1,1}, {3/2, 3/2}, to), (16) 

where p = R s /h, m = l/(4p 2 ), K(m) and E(m) are the 
complete elliptic integrals of first and second kind, re- 



spectively, and p F q ({ai 



}, ({&i, . . .,b q },m) is the 



generalized hypergeometric function. We consider ultra- 
thin film p > 1. Employing an appropriate expansions 
for the elliptic integrals and hypergeometric functions we 
find 



f 1 

C D 



2 



2 

7T p 

T28 



26 
225 



1 



15(V 



P » 1 (17) 
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FIG. 2: Components of the magnetic field of the tip I18H 
and 1191 at a distance d from the film: d/lo 0, 2, 5 and 10. 
Ha ~ /3Mo is the anisotropy field of the film and lo = 10 nm 
is the magnetic length. Solid and dashed lines correspond to 
H z and H r components, respectively. The radius of magnetic 
tip R = 10 Zq = 100 nm, the ratio of magnetization of the tip 
to that of the film is M t /M = 10. 



The expression 1)170 can be replaced by p 2 /2 - 



; p/128 



in the interval p > 1 with the accuracy of half percent. 

The spin - polarized scanning tunneling microscopy 
(SP-STM) and spectroscopy (SP-STS) offers a great po- 
tential to unravel the structure of a single Skyrmion and 
Skyrmion arrays on the nanometer scale. 



This value diverges logarithmically when the distance be- 
tween the tip and the magnetic film approaches zero: 



H r (R, d)«M t 



d< R 



7T R 3 /d 3 d > R 



In a special case when the distance between the film 
and the tip d = the components of the magnetic field 
are: 



H z (r) = -2irM t G[R-r], 



H r (r) = nM t (£) Vi (~, ~,2, (~) 



(21) 
(22) 



where a = sign(i? — r) and O is the Heaviside function. 
The expression (1221 logarithmically diverges at the dis- 
tances close to the tip's edge r = R 



H r {r) fa 2M t ln 



8R 



e 2 \r-R\ 



(23) 



The components of magnetic field of the tip versus r at 
different d are depicted in Fig. 

At large distances r>E the magnetic tip is equivalent 
to a magnetic charge Q = M t TrR 2 and asymptotically the 
magnetic field is Coulomb-like: 



H z (r,d) 
H r (r, d) 



ixM t R 2 d 
■KM t R 2 r 

( r 2 +d 2)3/2' 



(24) 
(25) 



To make tip stronger by increasing magnetic charge of 
the tip Q it is possible to increase the radius of the tip 
or its magnetization. The first method is more effective. 



4. CYLINDRICAL MAGNETIC TIP 

Consider a semi-infinite cylindrical magnetic tip. The 
stray magnetic field of the cylindrical tip is well known. 
Its z and r components are 

/>oo 

H z {r,d) - -2nM t R dq J^qR) J (qr) e~« d , (18) 
Jo 

poo 

H r (r,d) = 2nM t R / dq J^qR) J^qr) e~ qd , (19) 

where M t , R and d are the magnetization, the radius 
of cylinder and the distance between the tip and mag- 
netic film plane respectively. The radial component of 
the magnetic field of the tip H r (R, d) at the edge r = R 
can be expressed in terms of the hypergeometric function: 

H r (R,d)=nM t ^ 2 F 1 ^,3,-^pj, (20) 



5. STATIC EXCITATIONS IN THE PRESENCE 
OF THE TIP 

In order to study the relaxation dynamics of the 
Skyrmion created by the magnetic tip it is necessary first 
to analyze static solutions of the LLE for the magnetic 
film. A rough estimate of the minimal tip radius stabi- 
lizing the MBD at very small distance d can be extracted 
from the balance of energy: R$ = Sdw/ (27rMtMo), where 
£dw = 2y/af3MQ is the linear tension of the domain wall. 
For more accurate analysis the knowledge of the magne- 
tization and magnetic field distribution is necessary. 

Several cylindrically symmetric soliton-like solutions 
of the LLE have been found. They are the static 
Skyrmion solutionSii, the dynamic Skyrmion with the 
spin precession 8 and the slowly precessing magnetic bub- 
ble domain (MBD)2. 

The precessing solitons can live longer than solitons 
without precession, though they finally also collapse due 
to dissipation, in accordance with the Derick-Hobart 
theorem2ii22. The collapse does not proceed if the 
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Skyrmion configuration is stabilized by the higher order 
exchange interaction. Stabilization also occurs in ferro- 
magnets without center of inversion^. 



FIG. 3: Profile of 8(r) at d — and different values of the 
ratio of the tip magnetization to that of the film Mt/Mo. 
a) Profile of 9m bd for M t /M = 4,10. The radial compo- 
nent of magnetic field is shown on the same figure, b) A 
Skyrmion-like solution strongly distorted by the cusp of H r 
at the edge of the tip at Mt/Mo = 5,7,8. The energies of 
MBD □ and Skyrmion O versus Mt/Mo are shown in insert. 
Eo = 4-KahMo is the energy of the unperturbed Skyrmion. 



Let us rewrite the LLE in terms of spherical coordi- 
nates 9 and <fi, which determine the direction of magne- 
tization: 



— sin 0(9 = l 2 div (sin 2 9 (l + a 2 «(Vm) 2 )) - — sinfl sin(0 - <p) - — sin 2 9, (26) 

UJQ LOQ UJQ 

— sin 9 4> = Iq (- V 2 6> + (V0) 2 sin 9 cos 9) (1 + a 2 /i(Vm) 2 ) - Z 2 a 2 n div((Vm) 2 ) + 

sin 9 cos 9 + - — - sin 9 — - — - cos 9 cos(0 — tp) H 9, (27) 

ojq loo 



where (Vm) 2 = (V6») 2 + sin 2 9 (V0) 2 , ip is the azymuthal angle in the film plane; Iq — \J a,/ {5 is the magnetic length 
and tt>o = (3gMo is the ferromagnetic resonance frequency. In the static solution <fi = ip and 9(r) obeys an ordinary 
differential equation: 



d 2 9 ld9 _ . / 1 1 



-rr H ; smf cos„ , 

ar z r ar \ in r z 



« 2k f-4- fr(Vm) 2 ^ ) - (Vm) 2 ) - h z sin + h r cos 6 = 0, (28) 

yrar \ ar y H / 

where r = H z r /aMo. Equation l|28() should be solved with the following boundary conditions: 9(0) = tt, t9(oo) = 0. 

I 



One of the most powerful tools to solve the nonlinear 
ordinary differential equations is the fitting method. The 
heart of this method is the matching of a numerically 
derived solution to a known asymptotics at small and 
large distances from the center of nonlinear excitation by 
varying the shooting parameter. 



The key problem in this method is to find a proper 
trial function, with known asymptotics at r — * and 
r — > oo and fits it to the asymptotics of the real solution. 
We solve LLE for 9(r) (gSJ by the fitting method. We 



choose the trial function 9(r) in the following form: 



tan ■ 



0(f) R U) 



■ exp 



k 



(29) 



The trial function l|29|l has two advantages. First, it in- 
terpolates between two principal magnetic textures: the 
MBD and the Skyrmion by varying the fitting parameter 
Rv\ and second, it has a valid asymptotic at r — » 0. 
Being distorted by the radial component of the tip mag- 
netic field H ri the radial component of the magnetization 
acquires an unusual power asymptotics ~ 1 /r 2 . This be- 
havior is more pronounced when the tip is close to the 
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FIG. 4: The normalized energy of the MBD W/E J3TJ as a 
function of MBD radius Rb at different distances between the 
tip and film. Curves with □, O an d A correspond to d/R = 
0.1, 0.2, 0.5 respectively. The ratio of the radius of the tip to 
the magnetic lenth R/lo = 6 and fi = 5. 



magnetic film (d <C R). In this case the radial component 
of the tip magnetic field cannot be treated as a perturba- 
tion, since it causes a strong distortion of the Skyrmion. 
Small and large values of Rf correspond to the Skyrmion 
solution JSJ and the MBD, respectively. We match the 
numerical solutions of LLE with the known asymptotic 
at infinity. 

Let us now consider a special case when the distance 
between the surface of the and the tip d is zero. The angle 
0(r) at large distances decreases as 1/r 2 , similarly to h r . 
At small distances, the solution for 9 has the Skyrmion 
asymptotic. We have solved numerically equation 128|) 
for with the shooting method for different fj, — M t /Mo. 
For each value of /i we have found two different values of 
the fitting parameter Rf satisfying the matching condi- 
tion. We treat the solution with larger Rf as the MBD, 
whereas the solution with smaller Rf is interpreted as 
the Skyrmion-like solution. Profiles of 0(r) for different 
values of fi are presented in Fig. |21 The cusp of h r 
component of the magnetic field of the tip produces the 
noticeable distortion of Skyrmion solution in the region 
of the tip's edge. This distortion grows with fj, increasing 
as it is seen from Fig. [21 Strongly distorted Skyrmion- 
like excitation is the second solution of Landau-Lifshitz 
equation, but its energy is larger than that of the MBD. 
Thus, the magnetic tip nucleates MBD at static condi- 
tions. Fig. illustrates that the radius of the generated 
MBD equals roughly to the radius of the tip. This result 
can be obtained banalytically by the direct minimization 
of energy. 

In order to calculate the tip-film interaction en- 
ergy we substitute the expressions <|19ll f° r H g and 



H r and equation © into Using approximate 

relations: sech 2 ((r — R b )/l ) — 2l S(r — R b ) and 
sech((r — R b )/l ) — nl S(r — R b ), we find 



w z 



W z 
E 



-/ = -/ k(q)J 1 (qR)J 1 (qR b )e-' ld dq, (30) 



where k(q) = 2 n^R R b (1/q + n l /2) /(/? Z 2 ), E is the 
unperturbed Skyrmion energy: Eq — AitoiM^; and fi — 
M t /M - The total energy of the MBD w in units E is 
represented by the sum of 4 terms 



w 



we + w a +w z + w D , 



(31) 



where we, w a , wz and wd are the exchange, the 
anisotropy, Zeeman and demagnetization energies, re- 
spectively. The partial sum we + w a can be represented 
by an interpolation formula: 



w E 



lo Rb ' 



(32) 



valid with good precision in a broad interval of the ratio 
Rb/lo > 2. The first term in equation l|32|) is the linear 
tension of the domain wall, the second is the correction 
to it due nonzero topological charge of the MBD. The 
Zeeman energy is given by eq. (|30p . The demagnetiz- 
ing energy wd for an ultrathin magnetic film (h <C lo) 
slightly renormalizes the exchange energy and can be ne- 
glected in the first approximation. This approximation 
is correct provided the radius of the MBD is much larger 
than the thickness of the film. 

The energy of MBD versus its radius R b for different 
d is shown in Fig. 0] If d > d c = 0.385 R, the energy 
minimum at R b = R disappears; the tip cannot confine 
MBD anymore and the latter shrinks. 

Accepting R b = R, we find from Eq. I|31l) the range 
of parameters, in which the nucleation of the MBD is 
favorable: 



R 1 

r >- 



0.5 7rgg(p) 
U(P) ' 



where £ = 7T////3, 
9(p) 



3 3 
2F1 I -,-,3, 



m = 1 - 



2p_ 
ny/k 



[K(k) - E(k)} 



(33) 



(34) 



(35) 



p = d/R, k = (1 + p 2 /4)- x ; K(k) and E(k) are the 
complete elliptic integrals of the first and second kind, 
respectively. 

The total energy of the MBD has the following form 
w = C + 0.5/C- CC 2 f(p) cos^-0.57r£C<?(p), (36) 

where £ = R b /lo and ip is the global rotation angle. In 
the equlibrium state "0 = 0. 

At a fixed radius of the tip R, the bubble nucleates 
when the distance between the tip and the film becomes 
less than some critical value d r - 
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d/R 

FIG. 5: The phase diagram showing the range of parameters 
in which the nucleation of the MBD is energy favorable (region 
I), d is the distance between the tip and the surface of the film; 
R is the tip radius, fj, = 5, /3 = 100. 



6. DISCRETE SPIN MODEL WITH THE 
DIPOLAR INTERACTION 

In this section we study the dynamics of the bubble cre- 
ation and its transformation into a stable Skyrmion solv- 
ing numerically time-dependent Landau-Lifshitz equa- 
tion on a discrete lattice. We employ the model of a 
classical two-dimensional Heisenberg ferromagnet with 
the nearest, next nearest (NNN), and next to next near- 
est neighbor (NNNN) interaction and dipolar interaction. 
The corresponding Hamiltonian reads: 



H = -J^T S n -S n+a -0.5A^(S*) 

n, a ii 

jj ^ S n • S n ' ^ (S n • r nn ' ) (S n ' • r nn ' ) 



-J' S n ■ S n+ b + J" S„ ■ S n+ 2a 

n,b n,a 

]Ts n H, (37) 



interaction has the form F(q) = — (^/4)/(<7/2fcp), where 
v = pp m/(/T 2 h 3 ) is the density of states and 



i 



i 



2x 



-In 



1 



1 



(38) 



is the Lindhardt function, whose expansion into the Tai- 
lor series over powers of x is 



/(*) 



X 

y 



X 

15 



■2n 



(2n - 1) • (2n + 1) 



(39) 



Since both coefficients at x 2 and x , have the same sign, 
the ferromagnetic state provided by this interaction sat- 
isfies the condition of stability (see Section 0. 

In the continuous limit the Hamiltonian (|37|l takes a 
following form 



U = \ \ J(VSfd 2 x + Hj -:- H L , 



(40) 



where J = J + 2 J' — 4J", Hd is the demagnetization 
field energy. The perturbative part of the Hamiltonian 
Til is 



Hi 



J" a 4 



(V 2 S) 



A 



(l-S 2 )-hS 



d 2 x 



(41) 



where h = .g/isH is the scaled magnetic field. The 
Skyrmion is stable in the continuous model if both J and 
J" are positive. The discreteness changes the topology 
and make the continuous decay of Skyrmion possible. To 
study the discreteness effect in continuous model Haldane 
has proposed to remove one plaquette in the center of the 
Skyrmion 24 . An explicit calculation of the Skyrmion de- 
cay by this mechanism was performed in the workt The 
Skyrmion is metastable in a discrete lattice if in a con- 
figuration with only one reversed spin its rotation to the 
parallel orientation requires overcoming a potential bar- 
rier. This requirement is satisfied in a square lattice if 
J" > 0.15(J + 2J'). Though important for simulation on 
a discrete lattice, this effect does not play any role for 
real itinerant magnets, in which spin density is contin- 
uous. As we have shown earlier, the Skyrmion stability 
requires J" to be positive, whereas the sign of J' remains 
indefinite. 



where the summations run over all spin sites n, its near- 
est neighbors a, NNN along diagonals b and NNNN 
2a. We assume the ferromagnetic exchange J between 
NN, antiferromagnetic exchange J" between NNNN; the 
sign of NNN exchange J' is not fixed; A denotes easy 
axis anisotropy constant and D is the dipolar long range 
coupling constant. Both exchange interactions J' and 
J" mimic the real RKKY interaction in metallic fer- 
romagnets mediated by conduction electrons. In usual 
magnets the dimensionless ratio D/(Ja 3 ) is of the order 
10~ 4 — 10~ 3 . Fourier transform of the oscillatory RKKY 



7. FROM BUBBLE TO SKYRMION 

To study the dynamic transformation of the MBD into 
the Skyrmion we solved the time-dependent LLE by the 
Runge-Kutta method for a sample with free boundaries. 
The sample had a circular shape and contained 11308 
spins located at the sites of a square lattice. We choose 



time step 10 



for good convergence. For dynamical 



simulations we generated initially the MBD imitating a 
strong tip in a close vicinity of the magnetic film. The 
radius of the generated MBD is approximately equal to 
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FIG. 7: Magnetic bubble generated by numerical simula- 
tions. A small square distortion of axial symmetry is caused 
by symmetry of the lattice. This effect becomes pronounced 
during the dynamical transformation of the bubble into the 
Skyrmion. The ratio of the distance between the tip and film 
is d/R = 7/4. 



FIG. 6: Simulation of spin dynamics at different J' and J": 
crosses correspond to a stable Skyrmion, curve 1 corresponds 
to J" = 0.25 + 0.5 J', curve 2 to J" = 0.15 + 0.3 J'. Simula- 
tions are in a good agreement with the criterion of stability 
against long-wave spin perturbations. 



the radius of the tip (the difference between radii of the 
tip and the MBD is below 5 %). 

Once the bubble was generated, we started to move 
away the tip uniformly (d(t) = d + vt). The MBD 
gradually transformsed into a stable Skyrmion or col- 
lapsed depending on the exchange parameters J' and J" . 
The range of the exchange parameters corresponding to 
the Skyrmion stability is represented in Fig. |SJ The 
Skyrmion stability criterion is in good agreement with 
the numerical data. 

The spin structure of the MBD is shown in Fig. [7| 
The shape of the bubble is slightly distorted during the 
collapse. This distortion owes to the effect of the square 
lattice. 

We also studied the difference between continuous and 
discrete approximation by numerical simulations with 
different anisotropy constants A. The discrete models ap- 
proaches the continuous approximation when anisotropy 
constant A decreases and the domain wall width 1$ in- 
creases (See Table HJ. 

The process of transformation passes three stages. In 
the first stage the bubble domain is controlled by the 
removing tip until it reaches the distance d « 2R. During 
the second stage the MBD shrinks under the action of the 
linear tension. When its radius reaches the value R ~ lo, 
it can not be more described as the MBD and acquires 
the profile of a Skyrmion. A rough interpolation formula 



TABLE I: Normalized parameters of the soliton versus 
anisotropy constant A derived from numerical simulations. All 
parameters are normalized to ferromagnetic exchange inter- 
action J. Antiferromagnetic exchange interaction constant 
J" = 0.2 and J' = 0. 



Anisotropy 
constant, A 


0.001 


0.0025 


0.005 


0.01 


0.02 


0.03 


0.04 


Effective 
AFM 

exchange, J" e // 


0.19 


0.186 


0.181 


0.175 


0.166 


0.16 


0.155 


Skyrmion 
radius, R s 


4.28 


3.8 


3.22 


2.77 


2.36 


2.13 


1.97 


Domain 
wall-width, lo 


10.39 


8.78 


7.88 


7.36 


7.08 


6.79 


6.53 



for the transformation time T(i?) following from semi- 
intuitive consideration reads: 



v l v m T ujqT' 

In this equation the first term is the time during which 
the tip reaches the distance d = 2R from the film. At 
larger distance the interaction between the tip and the 
film is negligibly small. The second term corresponds to 
the shrinking of the MBD under the action of the Laplace 
force from the domain wall. The Laplace force is equal 

to Fl = cr/R, where a — 2"vJx/a is the linear tension. 
The local velocity of the domain wall is antiparallel to 
the radius and equal to T)Fl, where r\ = Tv m /(4X) is the 
magnetic viscosity. The time necessary to decrease 
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significantly the radius is equal to integral 



Tr(R) 



r dr 
o V ® ' 



This gives exactly the second term in equation (|42H . 
Finally, the third term is the relaxation time for the 
Skyrmion derived in£. This simple interpolation fits the 
simulation data with very good accuracy. 

A more sophisticated description of the transformation 
process can be obtained by substitution of the interpo- 
lation between the MBD and Skyrmion given by equa- 
tion l|29(l into the Hamiltonian Qand integrate it over all 
variables. The resulting expression is the Hamiltonian for 
slow variables R and ip. The dimensionless Lagrangian 
density I can be constructed given the Hamiltonian den- 
sity w by the following ansatz: 



I = 0.5 C 



2 W 

dr 



(43) 



where £ = R e /lo, r = u>ot and R e is the radius of topo- 
logical excitation. To incorporate the dissipation it is 
possible to use the Rayleigh dissipative function 7Z. In 
dimensionless variables it reads 1Z = ( (Ct + V ; ?)- We 
will not present the equations of motion here to avoid 
lengthy formulas with no essentially new results. They 
can be useful when a more detailed description of the 
process will be in order. 

For real magnetic film with Iq ps 10 nm, tip with radius 
R = 50 nm, damping constant r = 0.1 and tip's velocity 
v = lnm/ns transformation time is T ~ 125ns. The 
equilibrium size of the Skyrmion is 



Rq = a 



8J"\ 1/4 
3AL J 



1 -r- 3 nm, 



(44) 



where L = In [(3A/8J") 1/4 fo] and a is the lattice spac- 
ing. Skyrmion spin structure obtained by numerical sim- 
ulations of time-dependent Landau-Lifshitz Equation is 
represented in Fig. |SJ 

An effective magnetic field corresponding to the ex- 
change interaction J is Hj = J aj \xb- For Fe it is about 
89 T and for Co it is 112 T. For existing strong magnetic 
tips real value of created magnetostatic field is about 1 
T. Therefore we consider the ratio H/2irJ — 0.001. For 
iron film and tip with radius 50 nm minimum magnetic 
field of the tip H tip w 28 G. 



8. CONCLUSIONS AND DISCUSSION 

The main result of this work is the numerical evidence 
that the bubble domain of the 50 nm size initiated by a 
strong magnetic tip in a thin magnetic film shrinks into 
a very small, few nanometers metastable Skyrmion. For 
this purpose we employed the MC and the Runge-Kutta 
methods on the square lattice. The static Skyrmion and 
the magnetic bubble domain were found numerically as 
solutions of the static Landau-Lifshitz equation. 




FIG. 8: The Stable Skyrmion spin configuration generated 
from numerical simulations. The anisotropy constant is A = 
0.04, the exchange interaction constants are J = 1, J' = 0.6, 
J" = 0.5. 



The role of the magnetic tip is two-fold. Its magnetic 
field stabilizes the bubble domain, which otherwise is un- 
stable in a thin film. Being removed it does not confine 
the MBD anymore and the latter starts to shrink con- 
serving its topological charge. The Skyrmion is the fi- 
nal result of this shrinking. Thus, the position of the 
Skyrmion is determined with high accuracy by the tip. 
This is important for the experimental observation of the 
Skyrmion, since it is very small and its magnetic field 
is strongly localized as it was established by our calcu- 
lations. On the other hand, strong localization of the 
Skyrmion magnetic field makes possible to read the in- 
formation stored by Skyrmions if the distances between 
them and between them and empty spaces are more than 
3i? s , where R s is Skyrmion radius. Such a distance is 
sufficient to ensure that Skyrmions do not perturb the 
magnetization around. 

Limitations on exchange constants accepted in our 
work stem from the requirement of the Skyrmion stabil- 
ity. One of them is associated with the discrete character 
of the model and is not essential for a real itinerant mag- 
net with continuous distribution of magnetization. An- 
other one occurs in continuous model: the part of energy 
proportional to the fourth power of the gradient must be 
positive. Being imposed onto a lattice model it leads to 
an inequality, which provides sufficiently large antiferro- 
magnetic exchange between the NNNN. In real itinerant 
ferromagnets the natural RKKY interaction ensures the 
stability. 

Thus, the Skyrmion stability makes it is possible to 
create it by circular magnetic tip and to observe it by 
spin-polarized scanning tunnelling microscopy having a 
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nanometer resolution. It is not yet clear how the creation 
of the next Skyrmion by lateral scanning tip will influ- 
ence an already existing Skyrmion. When this problem 
will be solved, the Skyrmion arrays could be a promis- 
ing candidate for the magnetic memory of unprecedented 
high density above 10 TBit/in 2 . 
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